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Based on an axisymmetric galactic disk model, we estimate the equilibrium gas pressure P/k in 
the disk plane as a function of the galactocentric distance R for several galaxies (MW, M33, M51, 
M81, M100, M101, M106, and the SMC). For this purpose, we solve a self-consistent system of 
equations by taking into account the gas self-gravity and the presence of a dark pseudo-isothermal 
halo. We assume that the turbulent velocity dispersions of the atomic and molecular gases are fixed 
and that the velocity dispersion of the old stellar disk corresponds to its marginal stability (except 
for the Galaxy and the SMC). We also consider a model with a constant disk thickness. Of the listed 
galaxies, the SMC and M51 have the highest pressure at a given relative radius R/R25, while M81 
has the lowest pressure. The pressure dependence of the relative molecular gas fraction confirms the 
existence of a positive correlation between these quantities, but it is not so distinct as that obtained 
previously when the pressure was estimated very roughly [J], @|. This dependence breaks down for 
the inner regions of M81 and M 106, probably because the gas pressure has been underestimated 
in the bulge region. We discuss the possible effects of factors other than the pressure affecting the 
relative content of molecular gas in the galaxies under consideration. 
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1. INTRODUCTION 

Almost all of the active processes associated 
with star formation take place in a relatively nar- 
row gaseous layer embedded in the stellar disk. 
The gaseous layer is inhomogeneous in density and 
temperature and is more homogeneous in pressure, 
since a more tenuous medium is simultaneously 
hotter. For example, our Galaxy shows that the 
change in gaseous-layer thickness along the disk 
radius can be well reproduced by assuming a hy- 
drostatic equilibrium (see, e.g., @). The HI dis- 
tribution along the z coordinate indicates that the 
gas density at small z decreases approximately as 
exp(— z 2 /h1), where h z is the vertical scale height 
of the gas distribution, i.e., according to a law ex- 
pected for an isothermal gaseous layer inside a ho- 
mogeneous stellar disk, although there is a density 
excess at large z compared to this simple law [J, 

!■ 

Below, we will assume the galaxy to be axisym- 
metric and the gaseous disk to be in an equilibrium 
state in which its thickness is determined by tur- 
bulent gas velocities (although it is obvious that 
the equilibrium condition can be violated in the lo- 
cal regions associated with active processes in the 
disk). 

The mean equilibrium gas pressure at a given 
galactocentric distance R near the disk plane plays 
a very important (if not crucial) role in the tran- 
sition of the gas to molecular form @, 
H- Since the stars are formed in the molecu- 



lar gas layer, the star formation pattern and rate 
also depend on the conditions of the gas transi- 
tion from one phase to another, HI «-> H2. Having 
roughly estimated the turbulent equilibrium pres- 
sure of the gaseous layer in several disk galaxies, 
Blitz and Rosolowsky [l[, Q concluded that the 
relative molecular gas fraction increases almost lin- 
early with pressure. However, this important con- 
clusion needs to be tested. 

To calculate the pressure in galactic disks, the 
hydrostatic equilibrium and Poisson equations are 
commonly used and a number of simplifying as- 
sumptions are made. For an infinite disk with a 
vertical gas scale height much smaller than that 
for the stellar disk, when the contribution from 
the spheroidal components to the vertical potential 
gradient is disregarded, the pressure of the medium 
can be expressed by the formula (see, e.g., [l|) 

P = (2G) - 5 X gas V gas (p° s L + {\p 9 ast% (1) 

where G is the gravitational constant, £ gas is the 
total surface density of the gas, v gas is the velocity 
dispersion of the gas in z coordinate, and p s t a r and 
p gas are, respectively, the volume densities of the 
stellar and gaseous components in the disk mid- 
plane. For a selfgravitating isothermal stellar disk, 
the stellar surface density S s t ar = 2p s t a rh s tar, 
where the vertical scale height of the stellar disk is 
hstar = (v'itar/'ZnGpstar) ' 5 ■ Hence follows a sim- 
ple formula for the equilibrium turbulent pressure 
if the gas self-gravity is disregarded: 

P = 0.84(GS starf' '' 5 £ gas 7 q 5 ■ (2) 
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This expression is commonly used to calculate the 
pressure in the plane of galactic disks [H, 0], 



0. Note that the stellar disk thickness cannot 
be measured directly and it is generally taken 
rather arbitrarily. At the same time, v gas and 
h star are assumed to depend weakly on the galac- 
tocentric distance R and to change little from 
galaxy to galaxy. To a first approximation, the 
gas velocity dispersion is about the same in dif- 
ferent galaxies unless the regions of intense star 
formation or circumnuclear regions are consid- 
ered. In this case, the turbulent pressure is a 
function of only the stellar and gas surface den- 
sities, i.e., P{R) ~YP s f ar {R)Y, gas {R). In this ap- 
proach, not only the change in stellar disk thick- 
ness with R, but also the gravitational field of 
the gas and the dark halo is ignored. Blitz and 
Rosolowsky [l[ argue that this approach gives a 
pressure estimate with an accuracy of about 10% 
for T, s tar > 20 M Q /pc 2 . However, as we will show 
below, this error was underestimated significantly. 

In this paper, we estimate the equilibrium turbu- 
lent pressure of the interstellar medium that cor- 
responds to the mean gas density in the galactic 
plane at a given R for several disk galaxies. We 
obtain our estimates through a self-consistent so- 
lution of the equations that describe the vertical 
volume density distributions of the stellar, atomic, 
and molecular components of a disk that is as- 
sumed to be axisymmetric. We take into account 
the gas self-gravity, the possible change in stellar 
disk thickness with R, and the contribution from 
the dark halo to the gravitational potential of the 
galaxy. 



2. THE SAMPLE OF GALAXIES AND 
ADOPTED PARAMETERS 

To estimate the equilibrium gas pressure at var- 
ious R, we must specify the radial surface den- 
sity distributions of the main galactic components 
and the velocity dispersions Vi(R) for the stellar, 
atomic, and molecular disks. We chose several 
galaxies for which fairly complete data on the gas 
distribution, brightness, and rotational velocity are 
available in the literature: our Galaxy, M33, M51, 
M81, M100, M101, M106, and the Small Mag- 
ellanic Cloud (SMC). Basic characteristics of the 
galaxies are given in the table [J All of the esti- 
mates used were reduced to the given distances to 
the objects. 

For M81, M101, M106, and the SMC, the disk 
parameters were found by modeling the rotation 
curves (best-fit model); for the remaining galax- 
ies, the surface density distributions of the stellar 
disks either were taken from available publications 
or were recalculated from the radial brightness pro- 
file. The mass-to-light ratio calculated from stellar 
models based on disk color indices was used to 
pass from brightness to surface density. 

The estimation of the vertical stellar disk scale 



height (half-thickness) needed to determine the gas 
pressure presents a particular problem. For every 
galaxy, we used two models. In the first model, 
which we take as the preferred one, the stellar disk 
half-thickness for the sample galaxies (except for 
our Galaxy and the SMC) was estimated by as- 
suming that the velocity dispersion of the old stars 
constituting the bulk of the disk mass was minimal 
to ensure its dynamical stability. Strictly speak- 
ing, the half-thickness estimates obtained in this 
way are minimal. However, analysis of the obser- 
vational data shows that the velocity dispersion for 
most spiral galaxies is actually close to its critical 
value (see, e.g., [ll|, 0i ^ n tne second, sim- 

pler model, the stellar disk thickness was assumed 
to be constant along R. The pressure estimation 
procedure for this case is described in the next sec- 
tion. 

To determine the half-thickness of a marginally 
stable stellar disk, we initially found the epicyclic 
frequency from the model rotation curve: 



(3) 



where fi(i?) = V(R)/R is the angular velocity. 
For the gravitational stability of a collisionless, in- 
finitely thin homogeneous disk with respect to ax- 
isymmetric perturbations, the critical velocity dis- 
persion according to the Toomre criterion is 



Ccrit (R) 



3.36GS Btar (fl) 
ae(iZ) 



(4) 



In the general case, the critical radial velocity dis- 
persion is 



(v r )s 



QC C 



(5) 



The parameter Q in our paper is assumed to be 
constant along the radius and equal to 1.5. Sim- 
ulations of marginally stable collisionless galactic 
disks show that this value agrees well (to within 
~ 30%) with the numerical results in a large R in- 
terval — except for the central region where the 
bulge dominates and the outermost regions where 
Q can be twice as high (see, e.g., [HI). Note 
that an underestimation of Q and, hence, (v r ) s t ar 
means an overestimation of the gas density and 
pressure. 

As residual velocity measurements for old disk 
stars show, the ratio of the vertical and radial 
velocity dispersions for most spiral galaxies lies 
within the range 0.5-0.8; for earlier-type galaxies, 
this ratio is probably, on average, higher [lq |. [H ]. 
For all galaxies, we will use the approximate rela- 
tion: 



( v z)s 



0.5{v r 



(6) 



For the old disk of our Galaxy and the SMC 
disk, the vertical stellar velocity dispersions were 
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(10) (11) 


(12) 


(13) 


MW 


— 


20 


3.20 


640.9 


[3] 


[18] 


[18] 


220.0 


5.0 [17] 


— 


[19] 


M33 


0.7 


7 


1.18 


439.5 


[20] 


[21] 


[22] 


136.6 


7.0 [23] 


[24] 





M51 


8.4 


27.4 


4.38 




[27] 


[28] 


[28] 


120.0 


3.25 — 


[24] 





M81 


3.63 


11.55 


2.8 


1709.5 




[29] 


[28] 


88.0 


4.6 


[24] 





M100 


17.0 


18.3 


4.00 




[25] 


[25] 


[26] 


272.0 


4.7 


[24J 




M101 


7.48 


23.80 


4.6 


628.8 




[8] 


[8] 


236.0 


5.2 


[24] 




M106 


6.95 


10 


6.3 


933.7 




[28] 


[28] 


157.0 


8.0 


[24] 




SMC 


0.06 


3.47 


1.4 


168.0 




[3JJ 








[21] 


[32] 



TABLE I: Columns (2) and (3) give the assumed distance to the object in Mpc and its photometric radius 
D25/2 in kpc; (4) and (5) list stellar disk parameters: the radial disk scale length in kpc and the surface density 
corresponding to R = 0, £ s t ar (0), in units of Mq/pc 2 ; (9) and (10) list dark halo parameters: the asymptotic 
velocity in km/s and the core radius in kpc; all of the remaining columns give references to the papers containing 
the corresponding characteristics. 



taken from publications (see the table U for refer- 
ences). For the SMC, observations point to the 
velocity dispersion, (v z ) s t a r — 27.5 km/s, which 
changes little with galactocentric distance and is 
close to the velocity dispersion for the atomic gas 
(22 km/s). The latter value is atypically high for 
the galaxies, which is probably due to the inter- 
action between the Magellanic Clouds. Since the 
abundance of molecular hydrogen in this galaxy is 
low with respect to HI, the contribution from H2 
to the pressure was ignored. 

For all of the sample galaxies, except the SMC, 
we took constant values of (w z )hi — 9 km/s and 
(i , z)h 2 = 6 km/s. Other model parameters are 
the core radius R c and the central volume density 
Pdm(0) (or asymptotic velocity V as ) for a pseudo- 
isothermal galactic halo. For our Galaxy and M 33, 
these parameters were taken from published data 
(see the table); for the remaining spiral galaxies, 
these were found by modeling the rotation curve. 
The profile of a pseudo-isothermal halo is described 
by the formula 



Pdm (r) 



Pdm (0) 



1 + (r/R c y 
In this case, the circular velocity is 



(7) 



V C {R) = [47^^(0)^(1 - § arctan A)] - 5 . 

At large R, the velocity approaches its asymptotic 
value of 



V as = [4ttG Pdm (0)RX 5 , 



therefore 



Pum{r) 



Vas 1 

AttG (R 2 c +r 2 )' 



(8) 



(9) 



Following Narayan and Jog [3| , we took the core 
radius for our Galaxy to be R c — 5 kpc and the 



asymptotic circular velocity to be V as = 220 km/s 
(as estimated by Mera et al. [TtJ )■ The volume 
densities of the SMC components were calculated 
without the bulge and dark halo, whose masses in 
this galaxy are probably low [31] . In estimating the 
disk parameters, we used the SMC rotation curve 
obtained from HI data and corrected for the gas 
velocity dispersion (asymptotic drift), taken from 
Stanimirovic et al. [31( . 



ESTIMATION OF THE EQUILIBRIUM 
PRESSURE 



3.1. The System of Equations 

The thicknesses of the disk components, the 
midplane volume density of the gas, and the corre- 
sponding pressure were estimated from the atomic 
and molecular hydrogen surface densities, which 
were assumed to be known. To calculate the thick- 
ness of a three-component disk (stars, HI, and H2) 
in a general gravitational potential, we used the 
same method as that applied by Narayan and Jog 
[3] for our Galaxy. We took into account both the 
self-gravity of the individual components and the 
gravitational influence of the halo, which can be 
significant in the outer disk regions. The main sim- 
plification made in this case is that all disks are as- 
sumed to be coplanar and axisymmctric, while the 
HI and H2 layers are assumed to be isothermal, 
i.e., the gas velocity dispersions are constant with 
radius. Since the disk thickness is much smaller 
than the radial density scale length, we ignored 
the contribution from the radial density inhomo- 
geneity to the potential gradient in z. In a similar 
way, but using slightly different input parameters 
and approximations, the disk thicknesses were es- 
timated by Abramova and Zasov [331 ] for several 
galaxies common to ours. 
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FIG. 1: Radial variations of the volume densities of 
the three components in the galaxies of our sample. 
The densities are given in units of Mq/pc 3 and the 
Galactocentric distances are given in fractions of the 
photometric radii ifes- 



The basic hydrostatic equilibrium equation is 

{{vz)D dp, _^d<j> l d<j> d 
Pi dz f-f dz dz 1 1 

where p is the volume density, —d(f>/dz is the force 
per unit mass along the z axis, <fi is the correspond- 
ing gravitational potential, the index i pertains to 
one of the three disk components (stellar disk, HI, 
or H 2 ), and the index d pertains to the spherical 
halo. 

The Poisson equation for a thin axisymmetric 
disk is 



E 



if 



The combination of the last two equations leads 
to an expression for the volume density distribu- 
tion at a given galactocentric distance: 



Pi 



dz 2 ((«,)?> 



d 2 (j)d + J_ / dpi 
dz 2 pi \ dz 

(12) 

where the term in the first brackets corresponds to 
the potential of the three-component disk inside 
the halo. The influence of the galactic bulge on 
the disk thickness was disregarded, since the vol- 
ume density distribution, along with the velocity 
dispersions of the disk gas and stars in the bulge re- 
gion, are poorly known. Therefore, the central re- 
gions of the galaxies were excluded from our anal- 
ysis. 

The system of equations for the stellar and 
gaseous components was solved numerically by a 
fourth-order Runge - Kutta method. The two nec- 
essary boundary conditions in the z = midplane 
of the disk are 



(Po) 



and 



dpi 
dz 



0. 



(13) 



dz 2 



(11) 



The following normalization condition should be 
added to this: twice the integral of the volume den- 
sity over the z coordinate for each disk component 
should be equal to the surface (column) density 
Sj(i?), which was assumed to be known in the R 
interval under consideration. 

The volume density of each of the three compo- 
nents was found by the method of iterations. The 
first step is the simultaneous solution of Eq. (fT2|) 
for p s tar{z) with appropriate boundary conditions 
at zero gas densities. Subsequently, the problem is 
solved for pm( z ) using the stellar volume density 
obtained at the previous step. At the next step, 
the system is solved for pn 2 (z) with the known 
Pstariz)) and pm(z). Indeed, this is not enough, 
since the stellar disk was calculated without al- 
lowance for the influence of the gas, while molecu- 
lar hydrogen did not affect the atomic gas density 
estimate. Therefore, we solved the system of equa- 
tions successively four more times using the den- 
sities calculated in the previous iteration at each 
step. In this way, we obtained self-consistent so- 
lutions for each component in the general gravita- 
tional field and, hence, the volume density distri- 
butions along z at a given R. 



3.2. Pressure Variation Along the Galactic 
Radius 



The gas pressure was determined from the al- 
ready obtained solutions to the hydrostatic equilib- 
rium and Poisson equations. To within the errors 
of our iterative calculations (several percent), the 
equilibrium pressure of the interstellar medium ob- 
tained during the solution is equal to the dynamic 
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R, pc z, pc 

FIG. 2: Contour maps for the volume densities of the Galactic gaseous components in units of cm -3 : (a) for the 
atomic gas, (b) for the molecular gas, (c) the map of the molecular gas fraction Ph 2 /(phi+Ph 2 ), and (d) pressure 
variations in z coordinate at various Galactocentric distances; the heavy line represents the curve for the solar 
distance. 



pressure in the disk plane: 

P = Pdyn = PHI«HI + PH 2 Wh 2 . ( 14 ) 

where the contribution from elements heavier than 
hydrogen (the coefficient 1.38) was taken into ac- 
count in the gas densities. 

We compared the radial pressure variations cal- 
culated by our method and from the simplified de- 
pendence ((2]) used by Blitz and Rosolowsky [l| , @ ■ 
Since these authors took into account neither the 
radial variations in the thicknesses of the disk com- 
ponents nor the gas self-gravity, the simplified for- 
mula gives a systematic difference (underestimate) 
of the gas pressure compared to our results. The 
discrepancy increases with galactocentric distance 
and ranges from 30% in the inner galactic regions 
to more than 40% at large R, because the gas and 
the dark halo increase in importance in the outer 
disk regions. 

The resulting estimates of the radial volume den- 
sity variations for the stellar and gaseous com- 
ponents in the midplane are illustrated in Fig[TJ 
(a) stellar densities, (b) atomic gas densities, and 
(c) molecular gas densities. The SMC and MW 
have the lowest stellar disk volume densities, while 



M81 and M 106 have the highest ones. We empha- 
size once again that the above density estimates 
for all of the objects, except our Galaxy and the 
SMC, were obtained by assuming a marginal gravi- 
tational stability of the disk and, strictly speaking, 
give an upper limit for the density. 



3.3. Our Galaxy 

According to our calculations, the stellar vol- 
ume density near the solar orbit is 0.05 M Q /pc 3 , 
while the atomic and molecular gas densities are, 
respectively, 0.02 and 0.08 in the same units, cor- 
responding to a total density of all disk compo- 
nents w 0.15 M Q /pc 3 . The error in the stellar disk 
surface density, which is known with an accuracy 
no higher than 30%, introduces the largest uncer- 
tainty into the estimate. 

Our calculations allow the volume densities of 
the components to be estimated at any point with 
coordinates (R, z). Figure [2] shows contours of the 
volume densities pm (a) and pn 2 (b). In agreement 
with the observational data (see, e.g., [34]), the 
atomic gas in our model forms a layer that widens 
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FIG. 3: Radial equilibrium pressure distribution for all our galaxies. The solid and dashed lines represent, 
respectively, the distributions obtained for our main model and by assuming the stellar disk thickness to be 
constant. 



with galactocentric distance, although the model, 
naturally, cannot reproduce the observed warp of 
the gas layer at R > 12 kpc. The molecular layer 
is narrower than the atomic one due to its lower 
velocity dispersion. 

Figure [2J; presents a contour map that shows 
the distribution of the molecular gas fraction with 
respect to the total mass of the gas component, 
Ph 2 /(phi+Ph 2 )- We see from the figure that this 
parameter decreases rapidly with z at Galactocen- 
tric distances up to 7 kpc and considerably more 
slowly at larger distances. However, our model dis- 
regards the molecular-to-atomic gas phase transi- 
tion with increasing z. Previously, Imamura and 
Sofue [34[ discussed the possibility of this transi- 
tion attributable to ultraviolet radiation pressure 
and density variations. According to the model 
by these autors, which is based on measuring the 
vertical molecular gas density profile, the rela- 
tive molecular fraction ph 2 /(phi+Ph 2 ) decreases 
by half from its maximum value at a distance of 
~ 80 pc from the disk midplane, while according 
to our model, this decrease takes place at w 150 
pc. This confirms the possibility of a rapid, on a 
time scale of less than 10 7 , transition of the gas 
from molecular to atomic form and back. Numeri- 
cal simulations show [35| that these transitions can 
actually occur on a short time scale (several million 



years) in the presence of supersonic turbulence. 

Figure [2ji illustrates the turbulent gas pressure 
variations in z coordinate at distances from 4.5 to 
14.5 kpc in the disk plane at 2-kpc steps. The 
heavy line indicates the dependence for the solar 
distance (8.5 kpc). This scheme clearly illustrates 
a decrease in the vertical pressure gradient \dP/dz\ 
with Galactocentric distance. As a result, the gas 
pressure at small z is higher in the inner disk, but 
the reverse is true at large z and the pressure is 
higher in the outer Galactic regions. 



3.4. Radial Gas Pressure Profiles 

Parallel with the model of galaxies described 
above, where the stellar disk thickness changes 
along R and is determined by the condition for its 
dynamical stability, we used a simpler model with 
a constant disk thickness. Since the influence of 
the gravitational fields from the gas and the dark 
halo on the stellar disk thickness is ignored in this 
approach, the problem of estimating the densities 
and half-thicknesses of the stellar and gaseous com- 
ponents ceases to be completely self-consistent. 

As the law p s tar{z), we took the formula for an 
equilibrium gravitating isothermal disk 
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FIG. 4: Radial distributions of the volume densities for the three disk components of M33 (left) and their vertical 
scale heights (right). The dark-gray, light-gray, and lightest colors indicate H2, HI, and stars, respectively. The 
variations in (v z ) sta r by a factor of 1.25 relative to the value taken for a marginally stable disk (solid lines) 
correspond to the shaded regions. 



Pstar(z) = Pstar- (O)sedl 2 ( — — I , (15) 

V \' l z)sta7- J 

where the vertical disk scale height was assumed 
to be proportional to its radial scale length 

(h z ) star = = 0.2h R (16) 

for all galaxies, except the SMC, and 0.3hn for 
the SMC. The half-thicknesses, densities, and pres- 
sures of the gaseous components in the gravita- 
tional field of the stellar disk were determined at 
the same gas velocity dispersions as those taken in 
the first (main) model. 

The equilibrium pressure distributions (the de- 
pendence of P I k on the radial coordinate normal- 
ized to the optical radius of the galaxy, R/R25) 
are shown in Fig.[3]for the main model (solid lines) 
and for the model with a disk of constant thickness 
(dashed lines). When passing to the latter, the 
overall pattern of the dependences for most galax- 
ies, in general, changed little; the difference be- 
tween the pressure estimates does not exceed a fac- 
tor of one and a half (Alg(P/h) > 0.2), except for 
the inner region of M106, where Alg(P/h) w 0.4. 

In general, the range of gas pressures in the 
galaxies under consideration is almost two orders 
of magnitude. Of all the sample galaxies, the SMC 
has the highest pressure, although the central vol- 
ume density of the stellar component for this ob- 
ject is lowest (Fig. [T]). The high pressure is related 
both to the high gas content and to the high gas 
velocity dispersion in this galaxy. The lowest pres- 
sure is in the outer parts of our Galaxy. 

By varying the input model parameters, we in- 
vestigated the degree of their influence on the re- 
sulting distributions of the gas volume densities 



and pressures as well as the vertical scale heights 
of the stellar and gaseous components. 

The solution turned out to be most sensitive 
to the adopted velocity dispersion, primarily for 
the stellar disk (v z ) s t a r- For illustration, Fig. 0] 
shows the radial distributions of the volume den- 
sities for the three disk components and their ver- 
tical scale heights (the dark-gray, light-gray, and 
lightest colors correspond to H 2 , HI, and stars, re- 
spectively) calculated, as described above, in terms 
of the model of a marginally stable disk and after 
(v z ) s tar was increased and decreased by a factor 
of 1.25 (shaded regions) for the galaxy M33. The 
central volume density of the stellar component de- 
creases approximately by a factor of 1.5 with in- 
creasing (v z ) s tar and increases by the same factor 
with decreasing (v z ) sta r- However, the correspond- 
ing changes in gas density and pressure are smaller, 
no more than 20%. Obviously, the pressure de- 
creases with increasing {v z ) s t a r and increases with 
decreasing (v z ) sta r- 



3.5. Relationship Between the Pressure and 
Molecular Gas Content 

The pressure is assumed to play a dominant 
role in the atomic-to-molecular gas phase transi- 
tion (see the Introduction). Based on the shield- 
ing condition for molecular clouds in an UV radia- 
tion field, Elmegreen [|| found that the ratio of the 
molecular gas volume density to the total gas den- 
sity should be proportional to P 2 2 j~ 1 , where j is 
the UV radiation density. The simplified pressure 
estimates by Blitz and Rosolowsky ([l|, 0|) dis- 
cussed above led to the conclusion that the relative 
molecular gas fraction ij — Sh 2 /Shi oc p - 8 ^ - 9 for 
galaxies with various star formation rates. Sev- 
eral galaxies turned out to differ from others by a 
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FIG. 5: Left: the dependence between of the molecular gas fraction rj and the hydrostatic pressure calculated for a 
marginally stable disk; the dash-dotted line reflects the average dependence (without M 81 and M 106), the dotted 
line indicates the straight line corresponding to the correlation obtained by Blitz and Rosolowsky I2|. Right: the 
same dependence for the model with a constant disk thickness, the dashed line marks the inner regions of M 81 
and M 106 (see the text). 



higher rj at the same pressure. The authors sug- 
gested that this was related to the interaction of 
the galaxies with the ambient medium. One of 
these galaxies, M 100, belonging to the Virgo clus- 
ter, also enters into our paper. However, accord- 
ing to our calculations, it exhibits no anomalous 
behavior, which also put into question the men- 
tioned interpretation. 

Figure [5] shows the dependences rj(P/k) for the 
galaxies of our sample using the two pressure cal- 
culation methods described above: by assuming 
the stellar disk to be marginally stable (left) and 
assuming its constant thickness (right). Qualita- 
tively, the two models yield similar results. All 
of the galaxies, except M81 and M106, definitely 
fall on the general dependence: the amount of 
molecular gas compared to that of atomic one also 
increases along with equilibrium pressure of the 
medium. The dependence is best fitted by r\ ~ P k 
(dash-dotted line), where k = 1.34 ± 0.44 for the 
first model (without M 81 and M 106). For compar- 
ison, the dots in the figure mark the line that de- 
scribes the dependence taken from [l[, 77 ~ po.92 
The disagreement with our relationship can result 
not only from the rougher pressure estimates in the 
cited paper, but also from the inclusion of galac- 
tic regions with very high 77 > 10, which are vir- 
tually absent in our dependences. As was shown 
by Blitz and Rosolowsky [2j, a linear dependence 
more likely reflects the pressure estimation method 
than the physical relationship between the quanti- 
ties being compared for r\ > 2. The reason is that 
the pressure estimated from the simplified formula 
is assumed from the outset to be proportional to 
the atomic gas density and, at large i] > 2, to the 
molecular gas density, which varies over a much 
wider range than the HI density. 



The SMC is absent in the diagram: there is 
very little molecular gas in this galaxy and it is 
distributed highly nonuniformly. Because of the 
low heavy element abundance in the interstellar 
medium, the conversion factor that relates the CO 
line intensity to the number of H2 molecules on 
the line of sight for the SMC is probably much 
higher than that commonly assumed for spirals, 
which makes it difficult to determine the molecu- 
lar gas mass. It is variously estimated to be from 
several million M Q to 3 • 10 7 M Q @ , but it does not 
reach O.IMhi even in the latter case, which is an 
order of magnitude smaller than the expected one, 
given the high gas pressure (Fig. [3]). This discrep- 
ancy was also pointed out by Leroy et al. 0] . As in 
the case of M 33, which also lies below the average 
line in the rj{P/k) diagram, it would be natural to 
associate this peculiarity with the intensity of the 
UV radiation, which should play the most impor- 
tant role in these two galaxies due to the low dust 
content and active star formation. 

Note the unusual behavior of M 81 and M 106 in 
the rj(P/k) diagram. It stems from the fact that in 
the inner regions of these galaxies at distances of 
several kpc from their cores, the gas surface den- 
sity ceases to increase or even decreases toward the 
center (the corresponding segments of the curves 
for these galaxies in Fig. [5] are marked by dashes). 
For this reason, the resulting pressure related to 
both gaseous components remains low. However, 
since the H 2 density does not exhibit the same 
deep central "dip" as the HI density, the molec- 
ular gas fraction turns out to be high in this case. 
A similar peculiarity probably also takes place in 
the Andromeda galaxy (M31), which is not among 
the galaxies considered here. Observations show 
that the HI density in this galaxy also decreases 



toward the center in the inner disk, while the rel- 
ative molecular gas fraction increases; this effect 
probably cannot be explained by a change in con- 
version factor [36j. Note also that although the 
conversion factor can be lower in circumnuclear 
galactic regions due to the higher metallicity in the 
gas, the three mentioned galaxies do not stand out 
among the remaining ones by their relative oxygen 
abundance O/H [37|. The fact that the molecular 
gas content in the inner regions of M 33, where the 
gas exhibits an underabundance of heavy elements, 
nevertheless, follows the general dependence T](P) 
is indicative of the absence of a close correlation 
between the gas metallicity and the degree of its 
molecularization. Similar reasoning is given in Q 
for the galaxy IC 10 as well. 

Another factor that can affect the Eb abun- 
dance and can lead to a high 77 in the inner re- 
gions of M31, M81, and M106 is a low UV radia- 
tion density. At Galactocentric distances of several 
kpc, where rj is large, the intensity of the short- 
wavelength radiation in M31 and M81 is actu- 
ally low (see the GALEX images of the galaxies), 
but this explanation is probably invalid for M 106, 
where the inner region experiences a starburst. 

An underestimation of the gas pressure by a fac- 
tor of 3-10 due to the action of some factors ignored 
in our simple model can in principle be responsi- 
ble for the unexpectedly high relative molecular 
gas content in the inner regions of certain galax- 
ies. Since such galaxies have a large bulge (just 
as many of the galaxies with a ring-like HI distri- 
bution), it would be natural to associate the high 
relative H2 mass precisely with its presence. 

The bulge action on the gas is twofold. First, the 
bulge produces an additional force that compresses 
the gas in the disk plane (in our calculations, we 
took into account only the halo) . To a first approx- 
imation, the additional bulge-produced pressure is 



Pfr ~ PgasQzhga 



l<TgasV b 2 h gas /R 2 , (17) 



where g z — GMbh gas /R 3 is the z acceleration com- 
ponent attributable to the gravitational field of a 
bulge with mass Mf, and V& is the circular velocity 
component associated with the bulge. It is easy 
to verify that at h gas < 100 pc and V& < 200 
km/s, the pressure Pb/k will not exceed signifi- 
cantly 10 4 K/cm 3 . Therefore, including the bulge 
alleviates the problem only slightly, but does not 
solve it. 

The second possibility is the presence of a ther- 
mal hot plasma associated with the bulge in which 
a high gas temperature is supported by old stars 
(SNI explosions). The presence of a hot has in the 
bulges of spiral galaxies manifests itself mainly as 
soft X-ray emission. Although the observational 



data are so far rather scarce, a gas with a temper- 
ature (l-j-7) x 10 6 K was detected in the bulges and 
halos of several galaxies, including our Galaxy and 
M31 (see, e.g., [381 ] . [39(). Since the pressure inside 
the cool gaseous layer cannot be lower than the ex- 
ternal pressure, a hot gas will play a crucial role at 
a gas number density of 10~ 2 -r- 10~ 3 cm -3 in the 
bulge; this can explain the pressure underestima- 
tion in the models that disregard the influence of 
the ambient medium. Probably for the same rea- 
son, the molecular gas fraction in lenticular galax- 
ies within several kpc of the center is nevertheless 
high, despite the low gas content in the disk (see, 
e.g., 0). 

Another factor that is disregarded in axisym- 
metric disk models is the presence of spiral arms. 
The gas pressure inside of them is higher and the 
molecular gas concentration to the spiral arms is a 
well-established fact. The gas-compressing shock 
waves cause an increase in the fraction of molec- 
ular gas and give rise to giant molecular clouds. 
A tenuous molecular gas also exists between spiral 
arms, although it is more difficult to detect (for a 
discussion, see HU). It is not yet clear how strong 
the influence of the spiral pattern is on the degree 
of gas molecularization in the galaxy as a whole. 
However, in any case, it by no means always plays 
a significant role. Indeed, in contrast to the mor- 
phology of the spiral pattern, the integrated ratios 
Mh 2 /Mhi in spiral galaxies are almost indepen- 
dent of the morphological type or luminosity of 
the galaxies 42j . Spiral galaxies in which the bulk 
of the gas is in molecular form are encountered 
among the galaxies of all morphological types, ex- 
cept the latest ones (43|. Note that the gas-rich 
M 51 and M 100 with both an extended spiral pat- 
tern and a high relative molecular gas fraction r\ 
stand out among the galaxies considered in this 
paper by high fraction of molecular gas. However, 
the latter agrees well with the higher equilibrium 
gas pressure in these galaxies (see Fig. [SJ. 

Thus, we have confirmed the dependence of the 
molecular gas fraction on its mean equilibrium 
pressure at a given galactocentric distance, al- 
though it is not so distinct as the dependence ob- 
tained when P is estimated from very simplified 
formulas. However, such factors disregarded in the 
model as the external pressure on the disk and the 
intensity of the UV radiation can also play a sig- 
nificant role in some galaxies. 
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